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ABSTRACT 


Areas  of  improved  and  unimproved  soil  near  berthing  Pier  1 at  Treasure  Island, 
California,  were  investigated  by  the  Spectral-Analysis-of-Surface-Waves  (SASW)  test.  The 
upper  1 2 m of  sand  fill  beneath  the  approach  to  the  pier  had  been  densified  by  a vibrating  probe 
technique  in  1985.  The  area  of  improved  soil,  which  is  23  m wide  and  93  m long,  performed 
well  during  the  1989  Loma  Prieta  earthquake  (Mw  = 7.0),  while  sinkholes,  sand  boils  and 
cracks  formed  in  the  adjacent  unimproved  areas.  SASW  tests  were  conducted  on  a 240-m-long 
alignment  that  extended  across  the  area  of  improved  soil  using  a seismic  vibrator  as  the 
principal  source  with  receiver  spacings  from  7.6  m to  76  m.  Average  shear  wave  velocities 
determined  for  the  densified  and  undensified  sand  fill  below  the  water  table  were  192  m/s  and 
167  m/s,  respectively.  Two  simplified  analytical  procedures  based  on  shear  wave  velocity 
correctly  predict  no  liquefaction  for  the  densified  sand,  and  marginal  liquefaction  for  the 
undensified  sand.  Although  less  conservative,  liquefaction  assessment  procedures  based  on  the 
SPT  and  CPT  provide  similar  predictions. 


KEYWORDS:  building  technology;  SASW  test;  seismic  testing;  shear  wave  velocity;  soil 
improvement,  soil  liquefaction;  surface  waves 
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CHAPTER  1 


INTRODUCTION 


1.1  BACKGROUND 

Liquefaction  of  loosely  deposited  granular  soils  is  a major  cause  of  damage  in 
earthquakes.  Delineation  of  weak  soil  layers  and  prediction  of  their  liquefaction  potential  are 
key  inputs  in  the  engineering  design  of  new  and  retrofitted  structures.  This  information  is  also 
essential  for  reliable  estimation  of  economic  losses  during  future  earthquakes.  When  projects 
extend  for  great  distances,  such  as  lifelines  and  large  building  complexes,  cost-effective 
evaluations  of  extensive  areas  are  required.  Screening  techniques  based  on  geology,  hydrology, 
and  soil  conditions  show  promise  for  identifying  areas  requiring  more  rigorous  analyses. 
However,  even  these  areas  requiring  further  analyses  can  be  quite  large. 

One  promising  technique  for  spatially  evaluating  the  liquefaction  susceptibility  of 
granular  soils  is  the  Spectral- Analysis-of-Surface-Waves  (SASW)  test.  This  test  is  an  in  situ 
seismic  method  for  determining  small-strain  shear  wave  velocity,  Vs,  profiles  of  soil  deposits 
and  pavements  (Stokoe  and  Nazarian,  1985;  Stokoe  et  al.,  1988a;  Gucunski  and  Woods,  1991; 
Stokoe  et  al.,  1994).  The  SASW  test  does  not  require  boreholes,  and  has  the  advantage  of 
providing  broad  areal  coverage.  Testing  can  be  performed  at  sites  where  minimal  disturbance 
is  required  and  when  soils  are  difficult  to  sample.  The  use  of  Vs  as  an  index  of  liquefaction 
potential  is  justified  since  both  Vs  and  liquefaction  are  influenced  by  many  of  the  same  factors 
(e.g.,  void  ratio,  effective  confining  pressure,  stress  history,  and  geologic  age).  Thus,  the 
SASW  test  is  well  suited  for  profiling  large  areas  with  the  objective  of  developing  two-  and 
three-dimensional  Vs  images  of  the  subsurface. 

In  March  1996,  SASW  tests  were  conducted  across  an  area  of  densified  sand  at 
Treasure  Island,  California.  The  site,  called  the  Approach  to  Pier  site,  is  located  on  the  south- 
eastern corner  of  the  island,  as  shown  in  Fig.  1.1.  The  principal  objective  of  the  tests  was  to 
evaluate  the  ability  of  the  SASW  method  to  rapidly  delineate  stratigraphy  and  assess 
liquefaction  resistance  over  a significant  lateral  extent.  A partial  summary  of  the  test  results,  as 
well  as  liquefaction  assessment,  is  published  in  a conference  paper  by  Andrus  et  al.  (1998). 
This  report  presents  a complete  discussion  of  the  test  procedures,  analyses,  and  results. 


1 


Fig.  1.1  - San  Francisco  Bay  Showing  Locations  of  the  Approach  to  Pier  and  Fire  Station  Sites 
at  Treasure  Island. 


1.1.1  Treasure  Island 

Treasure  Island  is  a man-made  island  constructed  in  1936-37.  It  was  formed  by 
hydraulic  filling  behind  a perimeter  rock  dike.  The  perimeter  dike  served  to  contain  the 
hydraulic  fill,  and  was  raised  in  sections  over  the  previously  placed  fill.  Currently,  the  island  is 
owned  by  the  U.S.  Navy.  Plans  to  transfer  ownership  of  the  island  to  the  City  of  San  Francisco 
are  underway. 

Treasure  Island  was  classified  as  a level  I national  geotechnical  experimentation  site  at  a 
workshop  sponsored  by  the  National  Science  Foundation  and  Federal  Highway  Administration 
(Benoit  and  de  Alba,  1991).  Level  1 sites  are  those  sites  that  have  significant  national 
importance,  long-term  accessibility,  and  good  documentation.  Only  three  sites  were  classified 
as  level  I sites  at  the  workshop.  (The  other  two  sites  are  located  at  Texas  A&M  University.) 
The  purpose  of  national  geotechnical  experimentation  sites  is  to  provide  the  geotechnical 
engineering  community  a system  of  well-documented  sites  for  research  and  development. 

Much  of  the  work  to  date  at  Treasure  Island  centers  around  a full-scale  ground  response 
experiment  (de  Alba  et  al.,  1994;  de  Alba  and  Faris,  1996a).  Six  accelerometers  and  eight 
piezometers  are  operating  at  various  elevations  near  the  fire  station  (see  Fig.  1.1).  Inclinometer 
casings  are  in  place  at  the  fire  station  and  at  two  locations  along  the  perimeter  of  the  island, 
including  the  Approach  to  Pier  site. 
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1.2  PURPOSE 


Standard  procedures  for  conducting  the  SASW  test  have  not  been  developed.  The 
purposes  of  this  report  is  to  outline  a procedure  for  conducting  the  SASW  test  to  delineate  weak 
soil  layers  and  assess  liquefaction  susceptibility  at  sites  that  extend  for  great  distances.  It  is 
hoped  that  this  document  will  aid  the  engineer  in  planning  and  conducting  field  liquefaction 
investigations  using  the  SASW  test,  as  well  as  contribute  to  the  development  of  standard  test 
procedures. 


1.3  OVERVIEW  OF  REPORT 

Following  this  introduction,  in  Chapter  2,  the  SASW  field  testing  procedure  and  data 
analysis  are  reviewed,  and  methods  used  in  this  study  are  described.  The  field  investigations 
and  liquefaction  analysis  at  the  Approach  to  Pier  site  are  discussed  in  Chapter  3.  Chapter  4 
provides  a summary  of  conclusions  as  well  as  remarks  about  additional  needed  research. 
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CHAPTER  2 


SPECTRAL-ANAL YSIS-OF-SURFACE-WAVE  TEST 


2.1  INTRODUCTION 

The  Spectral- Analysis-of-Surf ace- Waves  (SASW)  test  is  an  in  situ  seismic  method  for 
determining  shear  wave  velocity  profiles  of  soil  deposits  and  pavements  (Stokoe  and  Nazarian, 
1985;  Stokoe  et  al.,  1988;  Gucunski  and  Woods,  1991;  Stokoe  et  ah,  1994).  The  SASW  test  is 
based  on  the  principle  that  high-frequency  (short-wavelength)  surface  waves  propagate  only  in 
the  near-surface  layers,  and  low-frequency  (long-wavelength)  surface  waves  propagate  through 
the  near-surface  layers  and  deeper  layers,  as  illustrated  in  Fig.  2.1.  Thus,  if  stiffness  varies  with 
depth,  surface  waves  of  different  frequencies  will  propagate  at  different  velocities. 
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Fig.  2.1  - Approximate  Distribution  of  Vertical  Particle  Motion  with  Depth  for  Two 
Surface  Waves  of  Different  Wavelengths  (after  Rix  and  Stokoe,  1989). 
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Surface  waves  have  been  used  to  study  soils  and  pavements  since  the  late  1940s.  As 
reviewed  in  Roesset  et  al.  (1991),  the  early  studies  generally  consisted  of  steady-state  harmonic 
vertical  excitation  using  heavy  shakers,  progressive  movement  of  receivers  to  determine  the 
resulting  surface  wave  wavelength,  and  analysis  by  simple  empirical  rules.  This  method  is 
known  as  the  steady-state  Rayleigh  wave  technique.  It  did  not  gain  wide  acceptance  because 
field  equipment  was  bulky,  test  procedures  were  cumbersome,  and  only  simple  sites  could  be 
correctly  determined  because  of  the  empirically-based  analysis.  During  the  last  two  decades, 
the  method  has  been  greatly  improved  with  the  development  of  personal  computers,  digital 
signal  analyzers,  simplified  test  procedures  to  obtain  velocity  measurements  at  several 
wavelengths  with  one  receiver  setup,  large  mobile  shakers,  and  accurate  computer  models 
which  provide  extensive  analysis  procedures.  These  new  techniques  are  collectively  known  as 
the  SASW  test. 

The  SASW  test  consists  of  three  phases:  (1)  field  testing,  (2)  data  reduction  and 
formation  of  the  experimental  surface  wave  dispersion  curve  (known  as  the  “signature”  of  the 
site),  and  (3)  inversion  of  the  experimental  dispersion  curve  to  obtain  the  shear  wave  velocity 
profile. 


2.2  FIELD  TEST  PROCEDURE 

The  basic  configuration  of  source,  receivers,  and  recording  equipment  in  the  SASW  test 
is  shown  in  Fig.  2.2.  Two  vertical  receivers  are  placed  on  the  ground  surface  a distance  D 
apart.  Surface  waves  are  generated  by  a source,  in  line  with  the  two  receivers,  at  a distance  D 
away  from  the  near  receiver.  The  source  can  be  some  type  of  vertical  impact,  a frequency 
sweep,  or  a random  vibration.  A FFT  signal  analyzer,  or  other  appropriate  device,  is  used  to 
record  the  two  receiver  signals,  and  then  transform  them  into  the  frequency  domain.  From  the 
two  frequency-domain  records,  the  coherence  and  the  phase  of  the  cross-power  spectrum  are 
computed.  These  four  records  are  saved  on  a floppy  diskette  or  a hard-drive  for  later  analysis. 


Computer  FFT  Analyzer 


Fig.  2.2  - General  SASW  Field  Testing  Configuration. 
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A typical  set  of  SASW  records  from  testing  performed  in  this  study  using  a seismic 
vibrator  is  presented  in  Fig.  2.3.  The  phase  of  the  cross-power  spectrum  (Fig.  2.3a)  represents 
the  phase  difference  between  signals  measured  at  the  receivers  as  a function  of  frequency.  It  is 
from  the  phase  of  the  cross-power  spectrum  that  travel  times  are  calculated  for  each  frequency, 
as  explained  in  Section  2.3.  The  coherence  function  (Fig.  2.3b),  ranging  between  0 and  1 for 
two  or  more  averages,  gives  a general  indication  of  the  quality  of  the  two  signals.  A coherence 
value  near  one  indicates  a very  high  signal-to-noise  ratio,  and  good  quality  data.  On  the  other 
hand,  a coherence  value  near  zero  indicates  a low  signal-to-noise  ratio,  and  likely  poor  quality 
data.  The  auto-power  spectra  (Figs.  2.3c  and  2.3d)  represent  the  distribution  of  relative  energy 
at  each  frequency,  and  are  a function  of  the  source  and  ground  characteristics.  A constant  auto- 
power spectrum  (Fig.  2.3c  between  frequencies  of  about  5 Hz  and  40  Hz)  indicates  constant 
relative  energy.  Variations  between  auto-power  spectra  for  Receivers  1 and  2 (Figs.  2.3c  and 
2.3d)  are  a result  of  ground  characteristics. 

According  to  Richart  et  al.  (1970),  about  two-thirds  of  the  energy  generated  by  a 
vertical  source  on  a uniform  half-space  is  transformed  into  surface  waves,  while  the  other  third 
is  transformed  into  body  waves  (compression  and  shear  waves).  In  addition,  body  waves  at  the 
surface  have  lower  amplitudes  than  body  waves  in  the  half-space.  Body  waves  also  exhibit 
geometrical  attenuation  in  a uniform  half-space  at  a rate  of  d'2  (where  d is  distance  away  from 
source),  whereas  surface  waves  geometrically  decay  at  a much  slower  rate  of  d'° 5 in  the  half- 
space. For  the  wavelengths  of  interest,  which  are  short  relative  to  the  long  travel  distances,  it  is 
reasonable  to  assume  that  the  body  wave  energy  is  generally  insignificant  compared  to  the 
surface  wave  energy.  Thus,  the  records  shown  in  Fig.  2.3  can  be  considered  to  be  mainly 
measures  of  the  surface  or  Rayleigh  wave  energy. 

Various  source-receiver  setups  and  instrumentation  requirements  for  SASW  testing  are 
reviewed  in  Gucunski  and  Woods  (1991).  The  preferred  source-receiver  setup  is  illustrated  in 
Fig.  2.4a.  In  this  setup,  called  the  common  receivers  midpoint  geometry , receivers  are  placed 
equal  distance  (D/2)  from  a fixed  centerline.  Testing  is  performed  with  the  source  on  both 
sides  of  the  receiver  array  to  compensate  for  effects  of  dipping  layers  and  any  phase  shifts  due 
to  receivers  and  instrumentation  (Nazarian  et  al.,  1983).  Testing  is  repeated  at  different  values 
of  D to  obtained  data  over  the  desired  range  of  wavelengths.  Unfortunately,  the  common 
receivers  midpoint  geometry  is  not  practical  for  automation  and  testing  of  large  areas. 

A source-receiver  setup  more  appropriate  for  automation  is  illustrated  in  Fig.  2.4b.  In 
this  setup,  called  the  common  source  geometry,  the  source  location  is  fixed  and  the  receiver 
centerline  location  is  varied.  Testing  is  performed  with  the  source  on  only  one  side  of  the 
receiver  array.  Hiltunen  and  Woods  (1990)  obtained  nearly  identical  results  from  tests  on 
pavements  using  the  two  geometries.  Nazarian  et  al.  (1995)  developed  a trailer-mounted  array 
of  receivers  and  source  in  the  common  source  geometry  for  testing  pavements.  In  their  trailer- 
mounted  system,  the  receivers  are  mechanically  lowered  onto  the  pavement  surface,  and  then 
raised  before  moving  to  the  next  test  location.  Moving  source-receiver  systems  for  generating 
continuous  shear  wave  velocity  profiles,  rather  than  profile  at  discrete  locations,  are  currently 
under  development  at  The  University  of  Texas  at  Austin. 
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a.  Phase  of  the  Cross  Power  Spectrum 


b.  Coherence  Function 
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Fig.  2.3  - Typical  Record  Set  Obtained  During  SASW  Testing  Using  a Seismic  Vibrator  as 
the  Source  at  the  Approach  to  Pier  Site;  Array  SA-6,  Receiver  Spacing  = 15.2  m. 
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Fig.  2.4  - Two  Source-Receiver  Setups  Used  in  SASW  Testing  (after  Heiesy,  1981). 
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2.3  EXPERIMENTAL  DISPERSION  CURVE  FORMATION 


The  cross-power  spectrum  is  defined  as  taking  the  Fourier  Transform  of  two  signals 
separately  and  multiplying  the  results  together.  From  the  phase  of  the  cross-power  spectrum,  a 
travel  time,  t,  between  receivers  is  obtained  for  each  frequency  by: 

t = O / (2  7t  f)  (2.1) 

where  O = phase  difference  in  radians,  and  f = frequency  in  Hz.  A Rayleigh  wave  phase 
velocity,  VR,  is  then  calculated  using: 

VR  = D/t  (2.2) 

The  corresponding  wavelength,  XR,  is  then  determined  from: 

= VR  / f (2.3) 

These  calculations  are  repeated  for  each  frequency  where  the  data  quality  is  good.  From  the 
results  at  all  receiver  spacings,  a plot  of  VR  versus  XR  is  assembled.  This  plot  is  called  the 
experimental  dispersion  curve.  This  curve  represents  the  “surface  wave  signature”  of  the  site 
and  is  the  culmination  of  the  field  measurements.  A sample  experimental  dispersion  curve  is 
shown  in  Fig.  2.5. 

2.4  SHEAR  WAVE  VELOCITY  PROFILE  DETERMINATION 

A shear  wave  velocity,  Vs,  profile  for  each  test  array  is  obtained  through  an  iterative 
process  of  matching  the  experimental  dispersion  curve  to  theoretical  dispersion  curves.  To 
begin  this  iterative  process,  called  forward  modeling , initial  properties  (shear  and  compression 
wave  velocities  and  densities)  and  layer  thicknesses  are  assumed.  Next,  a theoretical  response 
is  calculated  for  the  assumed  horizontally  layered  profile  using  three-dimensional  computer 
models  of  surface  wave  propagation  (Roesset  et  al.,  1991).  The  response  is  analyzed  to  obtain 
a theoretical  dispersion  curve.  The  theoretical  dispersion  curve  is  compared  with  the 
experimental  dispersion  curve.  Adjustments  are  made  to  the  assumed  elastic  properties  and 
layer  thicknesses  in  the  profile  until  satisfactory  agreement  between  the  theoretical  and 
experimental  dispersion  curves  is  obtained.  The  final  theoretical  dispersion  curve  for  the 
sample  experimental  curve  is  also  shown  in  Fig.  2.5.  It  is  assumed  that  the  elastic  properties, 
densities,  and  layer  thicknesses  for  the  final  theoretical  dispersion  curve  represent  the  actual 
profile  of  the  site.  Presented  in  Table  2.1  is  the  layered  model  used  to  determine  the  final 
theoretical  dispersion  shown  in  Fig.  2.5. 

Recent  work  by  Nazarian  et  al.  (1994)  and  Joh  (1996)  has  shown  that  the  forward 
modeling  process  can  be  automated,  at  least  for  simple  sites.  Nazarian  et  al.  report  that  the  time 
for  testing  and  determining  a Vs  profile  at  a site  like  the  one  shown  in  Fig.  2.4b  can  be 
accomplished  in  about  30  minutes. 
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Fig.  2.5  - Comparison  of  Experimental  and  Theorectical  Dispersion  Curves  from  SASW 
Testing  at  the  Treasure  Island  Approach  to  Pier  Site,  Array  SA-6.  (Note  that 
Array  SA-6  Lies  Within  the  Improved  Area.) 
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Table  2.1  - Layered  Model  Used  to  Determine  the  Final  Theoretical  Dispersion  Curve  for 
SASW  Array  SA-6  (see  Fig.  2.5)  at  the  Treasure  Island  Approach  to  Pier  Site. 
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2.5  SASW  PROCEDURE  FOR  THIS  STUDY 


Field  tests  were  performed  at  the  Treasure  Island  Approach  to  Pier  site  following  the 
general  procedure  outlined  in  the  preceding  paragraphs  using  a truck-mounted  seismic  vibrator 
(or  vibroseis)  weighing  18,000  kg  (40,000  lbs)  as  the  principal  source  and  velocity  transducers 
(or  geophones)  having  a natural  frequency  of  1 Hz  (Mark  Products1  Model  L-4)  as  the 
receivers.  The  vibroseis  (Mertz  Model  12),  shown  in  Fig.  2.6,  was  owned  and  operated  by 
Subsurface  Exploration  Company  of  Pasadena,  California.  During  testing  the  weight  of  the 
vibroseis  rested  on  a metal  plate  that  was  1.2  m wide  and  2.3  m long.  Conventional  vibroseis 
control  electronics  modified  at  The  University  of  Texas  at  Austin  were  used  in  conjunction 
with  a waveform  analyzer  (Hewlett-Packard  Model  3562A  Dynamic  Signal  Analyzer)  to 
control  the  vibroseis  output.  The  vibroseis  output  consisted  of  a sine  function  that  was  varied 
over  a range  of  frequencies.  At  higher  frequencies  the  sine  sweep  progressed  rapidly.  At  low 
frequencies  (less  than  about  5 Hz)  the  sine  sweep  progressed  rather  slowly.  Receiver  spacings 
used  were  as  long  as  76  m and  as  short  as  1.5  m.  Profiling  was  performed  by  marching  the 
source-receiver  array  along  an  alignment  using  one  or  two  spacings  (typically  15  m and  30  m). 
To  increase  productivity,  however,  hand-held  hammers  and  dropped  weights  weighing  about  10 
kg  (25  lbs)  were  generally  used  at  receiver  spacings  less  than  7.6  m.  Signals  from  each  receiver 
were  recorded  and  processed  using  the  analyzer.  A total  of  eight,  30-m-long  sections  were 
tested  using  various  receiver  spacings  at  the  site. 


Fig.  2.6  - Photograph  of  SASW  Testing  at  the  Treasure  Island  Approach  to  Pier  Site,  Array 
SA-2,  Using  a Vibroseis  as  the  Source  and  1 Hz  Geophones  as  the  Receivers. 


1 Trade  or  manufacture’s  names  appear  herein  because  they  are  essential  to  the  objectives  of  this  document.  The 
United  States  Government  does  not  endorse  products  or  manufacturers. 
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One  complicating  factor  at  the  site  was  the  asphalt  pavement  cover.  To  characterize  this 
stiff  layer,  additional  testing  was  performed  at  two  locations  (Arrays  SA-4b  and  SA-6b)  using 
accelerometers  having  natural  frequency  of  23  kHz  (Wilcoxon  Research  Model  728T)  and 
small  hand-held  hammers.  The  common  midpoint  geometry  was  used  for  these  tests.  Typical 
distances  between  accelerometers  were  0.15  m,  0.3  m,  0.8  m and  1.5  m.  Results  of  these 
additional  tests  were  used  to  establish  the  experimental  dispersion  curves  at  short  wavelengths 
(see  Fig.  2.5). 

To  obtain  shear  wave  velocity  profiles,  data  reduction  and  forward  modeling  programs 
developed  at  The  University  of  Texas  at  Austin  (Roesset  et  al.,  1991;  Joh,  1996)  were  used. 
Experimental  dispersion  curves  were  created  from  the  field  data  and  compacted  (or  averaged) 
using  the  computer  programs  called  WinSASW  and  Avg.Data7.  Theoretical  dispersion  curves 
were  created  and  matched  with  the  experimental  curve  for  each  receiver  spacing  using  a 
program  called  Fit7.  Fit7  allows  the  user  to  select  two-  or  three-dimensional  wave  propagation 
modeling.  The  three-dimensional  model  (Roesset  et  al.,  1991)  was  select  throughout  the 
forward  modeling  process.  The  assumed  properties  (primarily  Vs)  were  adjusted  until 
satisfactory  agreement  between  the  theoretical  and  experimental  dispersion  curves  was 
obtained.  Layer  thicknesses  were  based  on  available  borehole  and  cone  sounding  information. 
Where  no  layer  information  was  available,  such  as  the  asphalt  pavement  thickness  and  the  top 
of  bedrock,  layer  thicknesses  were  adjusted  to  obtain  the  best  comparison.  Agreement  between 
the  two  dispersion  curves  was  assessed  visually  and  by  a maximum  likelihood  method 
formulation.  Both  assessment  methods  were  used  throughout  the  forward  modeling  process. 
However,  the  maximum  likelihood  method  was  used  to  finalize  the  layered  model.  As  part  of 
the  maximum  likelihood  method  formulation  (Joh,  1996),  a root  mean  square  (rms)  error  was 
calculated  for  the  difference  between  the  compacted  experimental  dispersion  curve  and 
theoretical  dispersion  curve.  The  search  for  the  minimum  rms  error  involved  increasing  or 
decreasing  the  shear  wave  velocity  for  each  layer  in  the  assumed  profile  by  as  little  as  3 m/s. 
Presented  in  Appendix  A are  the  properties  and  layer  thicknesses  for  the  final  theoretical 
dispersion  curves.  Experimental  and  final  theoretical  dispersion  curves  for  each  SASW  test 
array  are  presented  in  Appendix  B. 


2.6  LIMITATIONS 

As  explained  by  Stokoe  and  Rix  (1987),  the  results  obtained  by  the  SASW  test  are 
considered  a more  averaged  measurement  than  are  results  from  other  in  situ  seismic  tests,  such 
as  the  crosshole  and  downhole  tests.  The  "averaging"  occurs  because  relatively  wide  receiver 
spacings  and  long  wavelengths  are  required  by  the  SASW  method  to  sense  material  properties 
at  depth,  compared  with  the  small  volume  of  material  sampled  in  the  crosshole  test.  A 
conceptual  view  of  this  difference  is  shown  in  Fig.  2.7.  The  two  tests  will  yield  similar  results 
if  the  material  properties  do  not  vary  laterally  at  the  site. 
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Source 


Fig.  2.7  - Conceptual  Comparison  of  Material  Sampled  During  S ASW  and  Crosshole 
Testing. 


The  depth  of  resolution  is  limited  by  the  receiver  spacing  and  the  ability  of  source  to 
generate  long  wavelengths.  Rix  and  Leipski  (1991)  suggest  that  the  maximum  wavelength  be 
at  least  1 to  2 times  the  desired  depth  of  resolution.  Results  obtained  during  this  study  indicate 
that  the  maximum  wavelength  obtained  is  about  2 to  3 times  the  receiver  spacing.  Thus,  the 
maximum  depth  of  resolution  is  roughly  equal  to  1 to  1.5  times  the  receiver  spacing.  This 
estimate  of  maximum  depth  of  resolution  is  reasonable  where  there  is  not  a large  velocity 
contrast  between  layers,  which  is  the  case  at  many  liquefaction  sites. 

Where  there  is  a large  velocity  contrast  between  layers,  such  as  soil  over  rock,  the 
maximum  wavelengths  should  be  3 to  10  times  the  depth  to  accurately  characterize  both  layers 
(Sohiel  Nazarian,  The  University  of  Texas  at  El  Paso,  personal  communication  to  Ronald  D. 
Andrus,  1998).  Thus,  receiver  spacings  may  need  to  be  as  much  as  5 times  the  depth.  This 
limitation  can  be  overcome  if  either  the  depth  to  rock  or  the  shear  wave  velocity  of  the  rock  is 
known.  One  approach  sometimes  used  to  estimate  the  shear  wave  velocity  of  rock  at  depth  is 
to  determine  the  compression  wave  velocity  of  the  rock  from  refraction  measurements  and 
assume  a Poisson’s  ratio.  While  not  necessary  for  liquefaction  evaluations  in  this  study, 
reported  field  shear  wave  velocity  measurements  made  at  other  locations  within  the  bedrock 
formation  were  assumed  in  order  to  establish  the  depth  to  bedrock  and  accurately  characterize 
the  overlying  stiff  soil  layer. 
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The  accuracy  of  the  SASW  inversion  process  depends  on  the  assumed  layering.  Rix 
and  Leipski  (1991)  found  that  if  the  assumed  layering  coincides  with  the  actual  layering  at  the 
site,  the  forward  modeling  process  will  generate  an  accurate  profile.  However,  if  an  assumed 
layer  spans  across  the  actual  layer  boundary,  a velocity  closer  to  the  average  for  the  two  layers 
will  be  determined.  This  limitation  can  be  avoided  by  using  layering  defined  by  other  field 
methods  or  by  using  layers  thin  enough  to  reproduce  the  variations  in  the  true  profile. 

The  uncertainty  in  the  value  of  shear  wave  velocity  determined  for  a given  layer  is  not 
easy  to  quantify,  and  current  analytical  procedures  are  dependent  upon  the  variables  of  the 
starting  model.  Rix  and  Leipski  (1991)  suggest  using  ratios  of  the  final  variances  to  the  initial 
variances  for  each  layer  as  an  indicator  of  uncertainty.  The  variances  are  a product  of  the 
inversion  process,  equal  to  the  square  of  the  standard  deviation.  They  found  that  at  a particular 
depth  the  ratios  (or  uncertainties)  are  greater  for  relatively  thin  layers  than  for  layers  that  are 
thick.  Joh  (1996)  suggests  evaluating  layer  sensitivity  by  examining  various  matrix  terms  (the 
row  component  of  the  model  parameter  spread  function  and  the  diagonal  term  of  the  resolution 
matrix),  and  adjusting  layer  thicknesses  to  obtain  layers  of  moderate  sensitivity.  The  root  mean 
square  error  between  the  experimental  and  theoretical  dispersion  data  is  then  used  in  the 
procedure  by  Joh  as  an  indicator  of  over  all  goodness  of  fit.  These  procedures  depend  on  the 
number  and  spacing  of  averaged  experimental  dispersion  curve  points. 

Two  limitations  of  using  Vs  measurements  to  evaluate  liquefaction  are  (Youd  et  al., 
1997;  Andrus  and  Stokoe,  1997):  (1)  Measurements  are  made  at  small  strains,  whereas 
liquefaction  is  a large  strain  phenomenon.  This  limitation  can  be  significant  for  cemented  soils, 
since  small-strain  measurements  are  highly  sensitive  to  weak  interparticle  bonding  which  is 
eliminated  at  large  strains.  It  can  also  be  significant  in  silty  soils  above  the  water  table  where 
negative  pore  water  pressures  can  increase  Vs.  (2)  No  samples  are  obtained  for  classification  of 
soils  and  identification  of  non-liquefiable  soft  clayey  soils.  Non-liquefiable  soils  by  the  so- 
called  Chinese  criteria  have  clay  contents  (particles  smaller  than  5 Jim)  greater  than  15%,  liquid 
limits  greater  than  35%,  and  moisture  contents  less  than  90%  of  the  liquid  limit  (Seed  and 
Idriss,  1982).  To  compensate  for  these  limitations,  a limited  number  of  borings  should  be 
drilled  and  samples  taken  to  identify  weakly  cemented  soils  that  might  be  liquefiable  but 
classed  as  non-liquefiable  by  Vs-criteria,  to  identify  unsaturated  silty  soils  that  might  have 
lower  values  of  Vs  should  the  water  table  rise,  and  to  identify  non-liquefiable  clay-rich  soils 
that  otherwise  might  be  classed  as  liquefiable. 

Finally,  large  ground  vibrations  generated  by  the  vibroseis  will  restrict  its  use  near 
structures.  In  this  study,  whenever  the  vibroseis  was  located  close  to  a structure,  ground 
vibrations  were  monitored  with  the  geophones  used  in  the  SASW  test.  The  vibration  level  was 
kept  below  about  2.5  mm/s  (0.1  in/s),  a factor  of  10  below  the  level  damaging  to  structures 
(Richart  et  al.,  1970).  For  structures  with  sensitive  equipment,  the  vibration  level  may  need  to 
be  kept  even  lower.  In  addition,  noise  produced  by  the  compressor  on  the  vibroseis  will  also 
restrict  its  use  in  residential  areas. 
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CHAPTER  3 


INVESTIGATION  AND  LIQUEFACTION  ANALYSIS 
AT  THE  APPROACH  TO  PIER  SITE 


3.1  INTRODUCTION 

The  Approach  to  Pier  site  is  located  on  the  south-eastern  comer  of  Treasure  Island  near 
Pier  1.  A photograph  of  the  area  is  shown  in  Fig.  3.1.  Lying  about  3 m above  mean  low  low 
water  (M.L.L.W.),  the  area  is  essentially  level  and  capped  by  a 127  mm  thick  layer  of  asphalt 
pavement.  At  the  water  front,  the  island  surface  slopes  into  the  water  at  an  angle  from  the 
horizontal  of  about  24  degrees.  This  sloping  surface  is  covered  with  rock  riprap. 

Because  of  concern  for  the  seismic  instability  of  the  water  front  slope,  the  fill  beneath 
the  approach  to  Pier  1 was  densified  to  a depth  of  about  1 2 m by  a vibrating  probe  technique  in 
1985.  As  shown  in  Fig.  3.2,  the  area  penetrated  by  the  vibrating  probe  was  23  m wide  and  97 
m long.  Following  the  1989  Loma  Prieta  earthquake  (moment  magnitude,  Mw  = 7.0),  no  signs 
of  ground  disturbance  were  observed  in  the  improved  area,  while  sinkholes,  sand  boils  and 
cracks  were  seen  in  the  adjacent  unimproved  areas  (Geomatrix  Consultants,  1990;  Mitchell  and 
Wentz,  1991).  In  March  1996,  Spectral-Analysis-of-Surface-Wave  (SASW)  tests  were 
performed  along  the  alignments  shown  in  Fig.  3.2.  General  subsurface  conditions  of  the 
improved  and  unimproved  areas  are  described  in  Section  3.2.  Liquefaction  effects  are 
discussed  in  Section  3.3.  In  Section  3.4,  results  of  SASW  testing  are  presented.  Simplified 
liquefaction  assessment  procedures  are  applied  in  Section  3.5. 


3.2  SUBSURFACE  CONDITIONS 

A generalized  two-dimensional  soil  profile  of  the  SASW  test  alignment  is  shown  in  Fig. 
3.3.  This  profile  was  developed  from  borings  and  cone  soundings  performed  by  other 
investigators  (Baker,  1985;  Geomatrix  Consultants,  1990;  Harding  Lawson  Associates,  1996). 
The  upper  12  m of  soil  is  sand  fill  initially  deposited  in  a loose  to  medium  dense  state  during 
hydraulic  filling.  Grain-size  distribution  curves  for  three  samples  taken  from  the  unimproved 
fill  are  shown  in  Fig.  3.4.  Samples  above  a depth  of  6 m contain  as  much  as  17%  fines  (silt  and 
clay).  Below  6 m,  samples  contain  1%  to  4%  fines.  The  fill  is  underlain  by  3 m of  native  silty 
clayey  sand  (Shoal  Sands)  followed  by  27  m of  soft  to  stiff  clay  with  interbedded  sand  layers 
which  thin  away  from  Yerba  Buena  Island  (see  Fig.  1.1).  The  clay  (Young  Bay  Mud)  is 
underlain  by  alternating  layers  of  very  stiff  sandy  clay  (Old  Bay  Clay)  and  dense  sand.  At  the 
time  of  SASW  testing,  the  water  surface  in  the  bay  was  about  2 m below  the  ground  surface. 


17 


Fig.  3.1  - Photograph  of  the  South-Eastern  Comer  of  Treasure  Island  Showing  the  Approach 
to  Pier  Site.  (Location  of  the  Vibroseis  is  Near  SASW  Test  Array  4b,  see  Fig.  3.2.) 
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Fig.  3.2  - Map  of  the  Approach  to  Pier  Site  Showing  Locations  of  Structures  (de  Alba 
and  Faris  1996a)  and  Tests. 
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Fig.  3.3  - Generalized  Soil  Profile  of  the  SASW  Test  Alignment. 
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Fig.  3.4  - Grain-size  Distribution  Curves  of  Three  Split-barrel  Samples  Taken  from  the 
Unimproved  Hydraulic  Fill  Sand. 
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Sandstone  and  shale  bedrock  occurs  at  a depth  of  87  m at  the  fire  station  (de  Alba  and 
Faris,  1996a).  It  is  assumed  that  the  bedrock  surface  slopes  upward  from  the  fire  station  to  the 
sandstone  rock  forming  Yerba  Buena  Island. 

3.2.1  Vibrating  Probe  Technique 

From  construction  drawings  by  Foundation  Contractor,  Inc.,  initial  vibrating  probes  and 
penetration  tests  were  conducted  at  the  northwest  comer  of  the  improved  area  to  determine  the 
optimal  probe  spacing.  Subsequent  production  probes  were  performed  to  produce  a final  1.90- 
m or  2.24-m  probe  spacing  in  a triangular  grid  pattern,  as  illustrated  in  Fig.  3.5. 

Mr.  Gerald  Manning  of  Foundation  Constructors,  Inc.  (personal  communication  to 
Ronald  D.  Andrus,  March  1998;  SI  conversions  are  shown  in  brackets)  described  the  vibrating 
probe  technique  as  follows: 

“The  tool  we  used  was  a KM6-6000E  Tomen  (Kencho)  Vibrator  which  had  a vibratory 
force  of  81  metric  tons. 

A circular  probe  was  clamped  to  the  Tomen  vibrator  measuring  [15  m]  fifty  feet  in  length. 
The  diameter  and  wall  thickness  was  [710  mm]  28”  x [16  mm]  5/8”  wall.  Elliptical 
shaped  windows  measuring  about  [150  mm]  six  inches  wide  by  [200  mm]  eight  inches 
high  were  cut  in  the  probe  and  spaced  about  [510  mm]  twenty  inches  center  to  center. 

The  specified  material  was  spread  over  the  area  to  be  improved  to  a depth  of  [0.9  m]  three 
feet.  The  material  was  pushed  up  against  the  probe  as  the  tube  was  vibrated  into  the 
ground.  The  natural  soil  offered  almost  no  resistance  to  the  penetration  of  the  probe.  To 
allow  the  material  to  enter  the  ‘windows’,  the  rate  of  advance  of  the  probe  was  slowed  to 
about  [2.5  m]  eight  feet  per  minute.  After  the  area  was  probed,  the  ground  was  back  at 
it’s  original  elevation. 

The  specified  material  was  a subgrade  (road  base)  mixture.  It  contained  about  three  to 
five  percent  fines.  It  is  the  writer’s  [Mr.  Manning]  opinion  that  we  should  have  achieved 
better  ground  improvement  had  the  material  been  free  of  all  fines.” 

Mr.  Manning  adds  that  the  metal-tube  probe  was  opened  at  the  bottom  to  permit  vibration  of 
the  specified  material  into  the  surrounding  sand  fill. 

Grain-size  distribution  curves  4,  5 and  6 shown  in  Fig.  3.6  are  for  three  split-barrel 
samples  taken  from  the  improved  soil.  Also  shown  are  the  curves  (1,2  and  3)  for  three  samples 
taken  from  unimproved  soil.  These  curves  suggest  that  the  specified  material  used  contained 
about  20%  to  40%  coarse  sand,  and  as  much  as  15%  gravel. 
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Fig.  3.5  - Construction  Drawing  of  Improved  Soil  Area  Showing  Locations  of  Vibrating  Probes  and  Penetration  Tests  Performed 
in  April  and  May  of  1985  (redrawn  from  construction  drawing  prepared  by  Foundation  Constructors,  Inc.,  1985; 
information  provided  by  Richard  Faris,  Naval  Facilities  Engineering  Command,  1996). 


Particle  Size,  mm 


Fig.  3.6  - Comparison  of  Grain-size  Distribution  Curves  of  Split-barrel  Samples  Taken  from  the 
Unimproved  and  Improved  Sand  Fill. 


The  criteria  for  improvement  was  based  on  the  Standard  Penetration  Test  (SPT)  and 
Cone  Penetration  Test  (CPT),  as  given  in  Table  3.1.  A multiplier  of  4.0  was  used  to  convert  the 
SPT  blow  count,  N,  to  equivalent  CPT  tip  resistance,  qc.  According  to  Mitchell  and  Wentz 
(1991;  1998),  the  SPTs  were  performed  at  points  equidistant  from  vibration  probe  locations 
and  depth  intervals  of  0.75  m.  The  average  of  three  consecutive  N values  was  to  be  no  less 
than  the  values  listed  in  Table  3.1.  The  CPTs  were  also  to  be  performed  at  points  equidistant 
from  vibration  probe  locations.  The  values  of  qc  were  to  be  no  less  than  the  values  listed  in 
Table  3.1,  except  where  the  friction  ratio  was  greater  than  2.0%.  Locations  of  boreholes  and 
cone  soundings  performed  during  construction  are  shown  in  Fig.  3.5. 
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Table  3.1  - Criteria  for  Soil  Improvement  (after  Mitchell  and  Wentz,  1998). 


Depth  below  ground 
surface,  m (ft) 

SPT  N value, 
blows  / 0.3  m 

CPT  qc  = 4N, 
MPa  (tons/ft2) 

1-5  ( 5) 

11 

4.2  ( 44) 

3.0  (10) 

15 

5.7  ( 60) 

4.6  (15) 

19 

7.3  ( 76) 

6.1  (20) 

22 

8.6  ( 90) 

7.6  (25) 

25 

9.6  (100) 

9.1  (30) 

27 

10.2  (105) 

10.7  (35) 

28 

10.9  (114) 

12.2  (40) 

30 

11.5  (120) 

Presented  in  Figs.  3.7  and  3.8  are  profiles  of  cone  tip  resistance  measurements 
performed  before  and  after  soil  improvement.  CPT  profiles  shown  in  Fig.  3.7  are  for 
measurements  performed  in  1985  at  the  time  of  soil  improvement.  CPT  profiles  shown  in  Fig. 
3.8  are  for  measurements  performed  in  1996.  Also  shown  in  these  figures  is  the  minimum  qc 
requirement  for  materials  with  friction  ratio  greater  than  2.0%.  Based  on  these  comparison,  the 
densification  criteria  appear  to  have  been  generally  met,  except  between  depths  of  3 m and  6 m 
where  the  sand  fill  contains  as  much  as  17%  plastic  fines. 


3.3  LIQUEFACTION  EFFECTS 

Geologists  from  Geomatrix  Consultants  and  the  United  States  Geological  Survey 
mapped  the  liquefaction  effects  which  formed  within  about  60  m of  the  slope  crest  during  the 
1989  Loma  Prieta  earthquake  (Geomatrix  Consultants,  1990;  Egan  and  Wang,  1991;  Power  et 
al.,  1998).  The  liquefaction  features  near  Pier  1 are  shown  in  Fig.  3.9.  It  is  interesting  to  note 
that  most  of  the  sinkholes  and  sand  boils  that  formed  lie  along  the  perimeter  of  the  improved 
area,  suggesting  that  maybe  the  ground  was  disturbed  and  loosened  by  the  vibrating  probes. 

Photographs  of  the  sinkhole  and  sand  boils  near  the  improved  area  are  presented  in  Figs. 
3.10  and  3.11.  Based  on  grain-size  analysis  by  Bennett  (1998),  the  sand  boil  shown  in  Fig. 
3.10b  contains  2%  gravel,  97%  sand,  and  1%  silt  and  clay.  The  median  grain  size,  D50,  and 
coefficient  of  uniformity,  Cu,  for  this  material  are  0.289  mm  and  1.9,  respectively.  These 
characteristics  are  very  similar  to  characteristics  exhibited  by  borehole  samples  2 and  3 take 
from  depths  of  6.8  m and  8.8  m shown  in  Fig.  3.6,  evidence  that  liquefaction  occurred  at  these 
depths. 
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Fig.  3.7  - Comparison  of  Criteria  for  Soil  Improvement  with  Cone  Tip  Resistance 

Measurements  Performed  Before  Improvement  (3  April  1985)  and  Shortly 
After  Improvement  (3  May  1985).  (information  provided  by  Richard 
Faris,  Naval  Facilities  Engineering  Command,  1996;  unpublished  cone 
penetration  test  data  by  V.  A.  Baker,  1985). 
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Fig.  3.8  - Comparison  of  Criteria  for  Soil  Improvement  with  Cone  Tip  Resistance 

Measurements  Performed  After  the  1989  Loma  Prieta  Earthquake  (15  May 
1996)  in  the  Improved  and  Unimproved  Areas  (Harding  Lawson  Associates, 
1996). 
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Fig.  3.9  - Map  of  the  Approach  to  Pier  Site  Showing  Locations  of  Structures  (de  Alba  and 
Faris,  1996a)  and  Liquefaction  Effects  Generated  by  the  1989  Loma  Prieta 
Earthquake  (modified  from  Geomatrix  Consultants,  1990;  Power  et  al.,  1998). 
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(a)  North  End  of  Sinkhole 


(b)  South  End  of  Sinkhole 

Fig.  3.10  - Photographs  of  Sinkhole  Near  Southern  Boundary  of  the  Area  of  Vibrating  Probes, 
as  Noted  in  Fig.  3.9.  (Photographs  by  Michael  J.  Bennett,  U.S.  Geological  Survey) 
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Fig.  3.1 1 - Photograph  of  Sand  Boil  Near  the  Area  of  Vibrating  Probes,  as  Noted  in  Fig.  3.9. 
(Photograph  by  Michael  J.  Bennett,  U.S.  Geological  Survey) 


Nearly  all  the  cracks  which  formed  in  the  asphalt  pavement  covering  the  unimproved 
areas  ran  parallel  to  the  water  front,  and  are  less  than  about  40  mm  wide  (see  Fig.  3.9). 
Summation  of  crack  widths  indicates  bayward  movement  of  the  perimeter  of  the  island  on  the 
order  of  80  mm  at  this  location.  This  amount  of  movement  is  relatively  small  compared  with 
the  several  meters  of  movement  that  can  occur  during  liquefaction.  Located  close  to  the 
perimeter  of  island,  sloping  ground  may  have  contributed  to  the  amount  of  lateral  movement. 
Based  on  these  findings,  liquefaction  in  the  unimproved  area  is  considered  marginal. 
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3.4  SASW  TEST  RESULTS 


SASW  tests  were  conducted  in  the  areas  of  improved  and  unimproved  soil  along  the 
240-m-long  test  alignment  shown  in  Fig.  3.2.  The  relative  locations  of  receiver  spacings  of  7.6 
m,  15.2  m,  and  30.5  m with  respect  to  the  area  of  vibrating  probes  are  shown  in  Fig.  3.12. 
Experimental  dispersion  data  obtained  for  these  29  test  setups  are  plotted  in  Figs.  3.13,  3.14  and 
3.15.  The  dispersion  data  for  test  arrays  in  the  improved  area  (solid  symbols)  are  distinctly 
separated  from  the  dispersion  data  for  test  arrays  in  the  unimproved  area  (open  symbols).  The 
dispersion  data  for  two  arrays  located  40%  to  50%  within  the  improved  area  (+  symbols)  lie 
between  the  open  and  solid  symbols,  as  shown  in  Figs.  3.14  and  3.15.  Values  of  VR  for  the 
improved  area  are  as  much  as  90  m/s  higher  than  values  of  VR  for  the  unimproved  area  at  a 
wavelength  of  3 m.  This  difference  in  VR-values  decreases  to  about  15  m/s  at  a wavelength  of 
30  m.  Between  wavelengths  of  5 m and  24  m,  the  average  difference  in  VR-values  is  31  m/s. 

SASW  test  array  4c  is  located  next  to  a sinkhole  formed  by  liquefaction  during  the  1989 
Loma  Prieta  earthquake,  as  shown  in  Fig.  3.2.  Values  of  VR  for  test  array  4c  are  among  the 
lowest  measured,  as  shown  in  Figs.  3.13  and  3.14.  This  observation  was  expected,  since  array 
4c  lies  closest  to  the  waterfront  slope  where  overburden  pressures  in  underlying  soils  are  lower. 
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Fig.  3.12  - Relative  Locations  of  SASW  Source  and  Receivers  With  Respect  to  the  Area  of 
Vibrating  Probes  for  Receiver  Spacings  of  7.6  m,  15.2  m,  and  30.5  m.  (Actual 
Location  of  the  Test  Alignment  is  Shown  in  Fig.  3.2.) 
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Fig.  3.13  - Experimental  Dispersion  Data  for  SASW  Tests  Conducted  at  the  Approach  to 
Pier  Site  with  Receiver  Spacing  of  7.6  m. 
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Fig.  3.14  - Experimental  Dispersion  Data  for  SASW  Tests  Conducted  at  the  Approach  to 
Pier  Site  with  Receiver  Spacing  of  15.2  m.  Dispersion  Data  for  Receiver 
Setup  Located  50%  with  the  Improved  Area  Indicated  by  + Symbol. 


32 


Surface  Wave  Velocity,  V , m/s 

R 

100  150  200  250 

1 


s 


& 

JS 

bi>  10 

c 

JD 

13 

> 

a 

£ 


100 


i I i r 


i i i r 


Test 

Array 


o 

□ 

> 

< 


2 

3 

4a 

4b 

5 

6 
7 


+ 8 


i — i — i — i — 

Receiver  Spacing  = 30.5  m 


J L 


Unimproved 

Area 


Improved  - 


40% 

Improved 

Area 


Fig.  3.15  - Experimental  Dispersion  Data  for  SASW  Tests  Conducted  at  the  Approach 
to  Pier  Site  with  Receiver  Spacing  of  30.5  m.  Dispersion  Data  for  Receiver 
Setup  Located  40%  with  the  Improved  Area  Indicated  by  + Symbol. 
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Shear  wave  velocity  profiles  for  SASW  tests  conducted  in  the  improved  and 
unimproved  areas  are  shown  in  Figs,  3.16  and  3.17,  respectively.  These  Vs-profiles  are 
determined  from  forward  modeling  using  experimental  dispersion  curves  composed  of  the 
dispersion  data  for  receiver  spacing  of  0.15  m to  61  m,  as  given  in  Appendix  B.  The  variability 
in  site  stiffness  indicated  by  the  Vs-profiles  is  similar  to  variability  in  site  stiffness  suggested  by 
the  dispersion  data  (see  Figs.  3.13,  3.14  and  3.15). 

Presented  in  Fig.  3.18  are  average  shear  wave  velocity  profiles  for  the  improved  and 
unimproved  areas.  Values  of  Vs  for  the  improved  area  are  about  94  m/s  higher  than  values  of 
Vs  for  the  unimproved  area  at  a depth  of  1 m (226  m/s  versus  132  m/s).  Between  depths  of  2 m 
and  13  m,  average  values  of  Vs  for  the  undensified  and  densified  fill  are  167  m/s  and  192  m/s, 
respectively.  In  other  words,  values  of  Vs  measured  in  saturated  soils  of  the  improved  area  are 
about  15%  higher  than  values  of  Vs  measured  in  the  saturated  soils  of  the  unimproved  area. 
This  increase  in  Vs  clearly  shows  the  improvement  to  the  site.  At  a depth  of  13  m,  the 
difference  between  Vs- values  is  about  7 m/s  (193  m/s  versus  186  m/s),  clearly  showing  little,  if 
any,  improvement.  The  trend  is  similar  to  the  measurements  of  VR  (see  Figs.  3.13,  3.14  and 
3.15).  Also,  the  trend  is  similar  to  the  values  of  N given  in  Figs.  3.16  and  3.17,  and  values  of 
cone  tip  resistance  given  in  Fig.  3.8.  A depth  of  13  m for  no  improvement  agrees  well  with  the 
reported  depth  of  densification  of  12  m. 

Assembling  the  dispersion  curves  and  Vs-profiles  presented  in  Figs.  3.13  through  3.17 
leads  to  the  two-dimensional  velocity  profiles  shown  in  Fig.  3.19.  Fig.  3.19a  was  constructed 
by  plotting  the  dispersion  curves  at  the  location  midway  between  receivers.  No  corrections 
were  made  to  the  dispersion  data  to  account  for  effects  of  receiver  spacing  or  to  make  the 
spatial  distribution  of  the  data  even.  Fig.  3.19b  was  constructed  by  plotting  the  Vs-profiles  at 
the  center  location  of  the  test  section  shown  in  Fig.  3.12  and  given  in  Appendix  A.  The  lateral 
limits  of  the  zone  of  vibrating  probes  shown  in  Fig.  3.19  are  from  construction  drawings  (see 
Fig.  3.5).  The  depth  limit  of  the  zone  of  vibrating  probes  is  12  m,  and  depth  is  roughly  1/3  to 
1/2  times  wavelength  based  on  an  empirical  rule  of  thumb.  Several  test  setups  near  the 
southern  end  of  the  improved  area  permit  good  resolution  of  the  boundary  separating  densified 
and  undensified  sands.  At  the  northern  end  of  the  improved  area,  however,  the  number  of  test 
setups  are  limited  and  the  agreement  between  the  velocity  profiles  and  the  lateral  limit  of 
vibrating  probes  is  rather  poor.  Nevertheless,  the  zone  of  densified  sand  is  clearly  identified  in 
both  velocity  profiles. 

Since  the  process  for  obtaining  dispersion  curves  is  not  computationally  intensive,  two- 
dimensional  profiles,  such  as  the  one  shown  in  Fig.  3.19a,  could  be  completed  during  field 
testing.  For  the  measurements  presented  in  this  report,  field  testing  was  completed  within  a 7 
hour  period.  This  time  could  be  reduced  once  a routine  is  established.  Thus,  similar  two- 
dimensional  profiles  with  lengths  of  500  m to  1000  m could  be  generated  in  a day. 

The  process  used  in  this  study  for  obtaining  a single  Vs-profile  was  computationally 
intensive,  often  requiring  more  than  8 hours  of  computer  time  to  complete. 
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Fig.  3.16  - Three  Shear  Wave  Velocity  Profiles  for  SASW  Tests  Conducted  in  the 
Improved  Area. 
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Fig.  3.17  - Seven  Shear  Wave  Velocity  Profiles  for  SASW  Tests  Conducted  in  the 
Unimproved  Area. 
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Fig.  3.18  - Comparison  of  Average  Shear  Wave  Velocities  Profiles  Determined  by  SASW  Testing  for  the  Improved  and 
Unimproved  Areas. 
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Fig.  3.19  - Two-dimensional  Velocity  Profiles  Showing  Zone  of  Improved  Soil. 

(Note  that  the  Amount  of  Data  Used  to  Construct  Profile  Shown  in 
Fig.  3.19a  Decreases  with  Depth,  see  Figs.  3.13-3.15;  and  the  Uncertainty 
of  Profile  Shown  in  Fig.  3.19b  Increases  with  Depth,  see  Figs.  3.16-3.17.) 
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3.5  LIQUEFACTION  ANALYSIS  USING  SIMPLIFIED  PROCEDURES 


As  reviewed  in  a recent  National  Center  for  Earthquake  Engineering  Research 
(NCEER)  workshop  paper  by  Andrus  and  Stokoe  (1997),  several  liquefaction  assessment 
procedures  based  on  Vs  have  been  proposed  during  the  past  decade.  Shown  in  Fig.  3.20  are 
five  curves  separating  liquefaction  and  no  liquefaction  proposed  by  various  researchers.  These 
curves  are  based  on  correlations  between  shear  wave  velocity  corrected  to  a reference 
overburden  stress,  VS1,  and  a cyclic  loading  parameter  called  cyclic  stress,  or  resistance,  ratio. 
A discussion  of  these  parameters  will  be  given  following  this  introduction.  The  curve  by 
Tokimatsu  et  al.  (1991)  was  developed  from  laboratory  cyclic  triaxial  test  results.  The  curves 
by  Robertson  et  al.  (1992),  Kayen  et  al.  (1992),  and  Lodge  (1994)  were  developed  from  limited 
field  performance  data.  The  curve  by  the  NCEER  workshop  was  developed  using  liquefaction 
and  non-liquefaction  case  histories  from  20  earthquakes  and  over  50  sites  in  soils  ranging  from 
sandy  gravel  with  cobbles  to  profiles  including  silty  clay  layers.  This  curve  is  considered  to  be 
a conservative  bound  for  uncemented,  Holocene-age  soils  with  fines  content  (particles  less  than 
75  pm)  of  5%  or  less,  since  average  values  of  VS1  for  the  critical  layer  were  used  in  its 
development. 

Shear  wave  velocity  and  penetration  measurements  from  the  Approach  to  Pier  site 
provide  an  important  opportunity  for  comparison  between  velocity-  and  penetration-based 
simplified  liquefaction  assessment  procedures. 

3.5.1  Liquefaction  Assessment  Based  on  VS1  and  CSR 

The  procedure  recommended  by  the  NCEER  workshop  (Youd  et  al.,  1997;  Andrus  and 
Stokoe,  1997)  follows  the  general  format  of  the  penetration-based  procedures,  where 
penetration  or  Vs  is  correlated  with  a parameter  called  cyclic  stress  ratio.  The  cyclic  stress 
ratio,  CSR,  at  a particular  depth  in  a level  soil  deposit  can  be  expressed  as  (Seed  and  Idriss, 
1971): 


CSR  = t Ja\  = 0.65  (a_/g)  (0/cr’v)  rd  (3.1) 

in  which  Tav  is  average  cyclic  shear  stress  generated  by  the  earthquake,  amax  is  maximum 
horizontal  ground  surface  acceleration,  g’v  is  initial  effective  vertical  (overburden)  stress,  Gv  is 
total  overburden  stress,  g is  acceleration  of  gravity,  and  rd  is  a shear  stress  reduction  factor  with 
a value  less  than  1.  Based  on  amax  of  0.16  g and  0.1 1 g recorded  in  the  x and  y directions  at  the 
fire  station  during  the  Loma  Prieta  earthquake  (Brady  and  Shakal,  1994),  an  average  value  of 
0.14  g is  assumed  for  the  analyses.  Overburden  stresses  are  calculated  using  density 
measurements  of  fill  materials  from  other  locations  on  Treasure  Island  with  similar  depth  and 
values  of  penetration  and  Vs.  The  density  of  the  unimproved  fill  is  about  1.76  Mg/m3  above  the 
water  table,  and  1.98  Mg/nr  below  the  water  table.  For  the  improved  fill,  the  density  is  about 
1.92  Mg/m3  above  the  water  table,  and  2.0-2.16  Mg/nr  below  the  water  table. 
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Fig.  3.20  - Comparison  of  Five  Liquefaction  Resistance  Curves  Based  on  Overburden 
Stress-Corrected  Shear  Wave  Velocity.  The  Curve  by  the  NCEER 
Workshop  (Youd  et  al.,  1997;  Andrus  and  Stokoe,  1997)  is  for  Uncemented 
Soils  of  Holocene  Age  with  Fines  Content  Less  than  or  Equal  to  5%. 


40 


The  shear  wave  velocity  is  corrected  with  respect  to  a reference  overburden  stress,  Pa, 
by  (Robertson  et  al.,  1992): 

Vs,  = Vs(Pa/a’v)0-25  (3.2) 

where  Pa  is  typically  100  kPa  or  approximately  atmospheric  pressure,  and  a’v  is  the  effective 
overburden  stress  in  kPa.  Equation  3.2  ignores  the  effects  of  horizontal  stress  on  Vs. 

At  Treasure  Island,  liquefaction  in  the  unimproved  soils  most  likely  occurred  where 
values  of  Vs,  are  least,  and  where  values  of  cyclic  stress  ratio  are  greatest.  These  conditions 
occur  between  the  depths  of  6 m and  12  m. 

Resistance  to  liquefaction  can  be  approximated  by  (Andrus  and  Stokoe,  1997;  modified 
from  Dobry,  1996): 


CRR  = a(Vsl/100)2  + b[l/(Vslc-  VS1)  - l/Vslc]  (3.3) 

where  CRR  is  the  cyclic  resistance  ratio,  Vsic  is  the  limiting  value  of  VS1  which  separates 
contractive  and  dilative  behavior,  and  "a"  and  "b"  are  curve  fitting  parameters.  The  assumption 
of  a limiting  value  of  VS1  is  equivalent  to  the  assumption  commonly  made  in  the  SPT-based 
procedure  where  liquefaction  is  considered  not  possible  above  a N-value  of  30  (see  Fig.  3.23). 
The  relationship  developed  as  part  of  the  NCEER  workshop  for  magnitude  7.5  earthquakes  and 
clean  sands  is  defined  by  Vsic  = 220  m/s,  a = 0.03,  and  b = 0.9.  For  magnitude  7 earthquakes, 
Eq.  3.3  is  multiplied  by  a scaling  factor  of  about  1.25. 

Using  Eq.  3.3  with  the  above  mentioned  best-fit  values,  the  curve  separating 
liquefaction  and  no  liquefaction  for  magnitude  7 earthquakes  is  shown  in  Fig.  3.21.  Also 
plotted  in  Fig.  3.21  are  values  of  VS1  and  cyclic  stress  ratio  for  the  critical  layer.  The  data  point 
for  the  improved  area  correctly  lies  in  the  region  of  no  liquefaction  defined  by  the  curve.  For 
the  unimproved  area,  the  data  point  lies  on  the  curve.  Thus,  marginal  liquefaction  is  correctly 
predicted  for  the  unimproved  area. 
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Fig.  3.21  - Comparison  of  the  Liquefaction  Assessment  Chart  Based  on  VS1  and  CRR 

Recommended  by  NCEER  Workshop  (Youd  et  al.,  1997;  Andrus  and  Stokoe,  1997) 
for  Clean  Soils  with  Results  from  the  Approach  to  Pier  Site  Between  Depths  of  6 m 
and  1 2 m. 


3.5.2  Liquefaction  Assessment  Based  on  Vs  and  amax 

Another  method  relating  liquefaction  potential  and  Vs  has  evolved  from  the  strain 
approach  by  Dobry  et  al.  (1982)  and  analytical  studies  by  Stokoe  et  al.  (1988b).  In  this  method 
liquefaction  resistance  is  related  to  Vs  and  amax  directly.  By  combining  Eqs.  3.1,  3.2  and  3.3,  a 
relationship  for  magnitude  7.5  earthquakes  based  on  Vs  and  amax  can  be  obtained  in  the  form  of 
(Andrus  and  Stokoe,  1997): 

Kjt  = f,{a(f2Vs/100)2  + b[l/(Vslc  - f2Vs)  - 1/Vslc] } (3.4) 

where  fj  = a’v/(0.65  av  rd)  and  f2  = (Pa/  a’v)0  25.  Assuming  (1)  the  water  table  is  located  midway 
between  the  ground  surface  and  the  center  of  the  most  vulnerable  layer  and  (2)  the  density  of 
soil  is  1.76  Mg/m3  above  the  water  table  and  1.92  Mg/m3  below  the  water  table,  then  fj  and  f2 
can  be  approximated  by: 


f,  = l.l/rd 

(3.5) 

f2  - (7.3/z)0'25 

(3.6) 

where  z is  the  depth  to  the  center  of  the  most  vulnerable  layer  in  meters.  Eqs.  3.4,  3.5,  and  3.6 
provide  a simple  relationship  between  Vs  and  amax  that  depends  on  depth. 
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Fig.  3.22  - Comparison  of  Liquefaction  Assessment  Chart  Based  on  Vs  and  a^  (Andrus  and 
Stokoe,  1997)  for  Clean  Soils  with  Results  from  the  Approach  to  Pier  Site  Between 
Depths  of  6 m and  12  m. 


The  curve  for  magnitude  7 earthquakes  and  depth  of  9 m is  shown  in  Fig.  3.22. 
Liquefaction  behavior  predicted  by  this  method  is  similar  to  the  method  based  on  VS1  and  CSR 
(see  Fig.  3.21). 

3.5.3  Liquefaction  Assessment  Based  on  (Nj)60  and  CSR 

The  most  widely  used  simplified  procedure  for  assessing  liquefaction  resistance  is  the 
procedure  developed  by  the  late  Prof.  H.B.  Seed  and  his  colleagues  (1971,  1982,  1983,  and 
1985)  based  on  modified  SPT  blow  count.  The  NCEER  sponsored  workshop  (Youd  et  al., 
1997)  reviewed  the  procedure  and  recommended  some  revisions.  This  updated  procedure  is 
applied  here  using  the  measured  standard  penetration  resistances  given  in  Figs.  3.16  and  3.17. 

The  measured  standard  penetration  resistance,  Nm,  is  modified  to  an  equivalent  clean- 
sand  standard  penetration  resistance,  (Nj)60,  by  the  following  equation: 

(N,)6„  = NmCNCECBCRCs  (3.7) 
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where  CN  is  a correction  factor  for  overburden  pressure,  CE  is  a correction  factor  for  hammer 
energy  ratio,  CB  is  a correction  factor  for  borehole  diameter,  CR  is  a correction  factor  for  rod 
length,  and  Cs  is  a correction  factor  for  samplers  with  or  without  liners.  The  overburden 
correction  factor  is  commonly  calculated  by  (Liao  and  Whitman,  1986): 

CN  = (Pa/a’v)05  (3.8) 

where  Pa  is  atmospheric  pressure  (about  100  kPa),  and  a\  is  the  effective  overburden  stress  in 
kPa.  The  Nm-values  given  in  Figs.  3.16  and  3.17  were  measured  by  Harding  Lawson 
Associates  (1996)  in  a 150-mm  borehole  using  a standard  sampler  with  a 35  mm  inside 
diameter  not  needing  liners  and  a safety  hammer.  Typical  correction  factors  for  this  equipment 
are:  CE  = 1.0,  Cg  = 1.05,  and  Cs  = 1.0.  For  the  rod  length  correction,  NCEER  workshop 
recommended  values  of  CR  are  0.95  for  lengths  of  6 m to  10  m,  and  1.0  for  lengths  of  10  m to 
20  m (Youd  et  al.,  1997). 

Plotted  in  Fig.  3.23  are  (N^-values  and  cyclic  stress  ratios  for  the  critical  layer 
between  depths  of  6 m and  12  m.  Also  plotted  is  the  boundary  separating  liquefaction  and  no 
liquefaction  recommended  for  magnitude  7 earthquakes  by  the  NCEER  workshop  (Youd  et  al., 
1997),  assuming  a magnitude  scaling  factor  of  1.25.  The  data  for  the  improved  area  correctly 
lie  in  the  region  of  no  liquefaction.  At  an  average  modified  blow  count  of  54  and  a minimum 
modified  blow  count  of  4 1 , these  data  points  lie  well  above  the  critical  value  of  30  suggested  by 
Seed  et  al.  (1985),  as  shown  in  Fig.  3.23. 


Modified  SPT  Resistance,  (Nj)^,  blows  per  0.3  m 

Fig.  3.23  - Comparison  of  Liquefaction  Assessment  Chart  Based  on  (Nj)60  and  CSR 

Recommended  by  the  NCEER  Workshop  (Youd  et  al.,  1997;  modified  from  Seed 
et  al.,  1985)  for  Clean  Sands  with  Results  from  the  Approach  to  Pier  Site,  Borings 
IA1  and  IA2  Between  Depths  of  6 m and  12  m. 
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For  the  unimproved  area,  the  average  and  minimum  values  of  (Nj)60  also  lie  in  the 
region  of  predicted  no  liquefaction.  With  boring  BH-IA1  located  close  to  a sinkhole  (see  Figs. 
3.2  and  3.9),  a prediction  of  no  liquefaction  appears  to  be  inconsistent  with  field  behavior.  One 
possible  explanation  for  this  inconsistency  is  that  liquefaction  was  localized  or  marginal  at  this 
site.  In  the  original  chart  by  Seed  et  al.  (1985)  and  in  a later  chart  by  Youd  and  Noble  (1997), 
there  are  a few  reported  liquefaction  and  marginal  liquefaction  case  histories  that  plot  in  the 
region  of  no  liquefaction.  Another  explanation  is  that  the  ground  shaking  caused  the  soil  to 
settle  and  densify,  resulting  in  higher  post-earthquake  penetration  and  velocity  measurements. 

3.5.4  Liquefaction  Assessment  Based  on  qclN  and  CSR 

In  recent  years,  the  Cone  Penetration  Test  (CPT)  has  become  popular  due  to  its  greater 
repeatability  and  nearly  continuous  profile.  The  abundance  of  field  performance  data  has  lead 
to  the  development  of  various  liquefaction  assessment  procedures.  These  procedures  were 
reviewed  during  a recent  NCEER  workshop  (Youd  et  al.,  1997;  Robertson  and  Wride,  1997; 
Olsen,  1997).  The  procedure  recommended  by  the  workshop  is  applied  here  using  the  CPT 
profiles  listed  in  Figs.  3.7  and  3.8. 

The  measured  cone  penetration  tip  resistance,  qc,  is  corrected  and  normalized  by  the 
following  equation: 


4c IN  — Cq  (4/^a) 


(3.9) 


where 


CQ  = (PJ&J  • (3.10) 

CQ  is  a correction  factor  for  overburden  pressure,  Pa  is  atmospheric  pressure  (about  100  kPa), 
and  n is  an  exponent  dependent  on  grain  size  characteristics  of  the  soil.  A typical  value  of  n is 
0.5  for  clean  sands. 

Plotted  in  Fig.  3.24  are  average  qclN-values  and  cyclic  stress  ratios  for  the  critical  layer 
between  depths  of  6 m and  12  m.  Also  plotted  is  the  boundary  separating  liquefaction  and  no 
liquefaction  recommended  by  the  NCEER  workshop  (Youd  et  al.,  1997)  for  magnitude  7 
earthquakes,  assuming  a magnitude  scaling  factor  of  1.25.  The  data  for  the  improved  area 
correctly  lie  in  the  region  of  no  liquefaction. 

For  the  unimproved  area,  the  average  values  of  qclN  also  lie  in  the  region  of  predicted  no 
liquefaction,  with  the  exception  of  one  lying  on  the  curve.  Thus,  marginal  to  no  liquefaction  is 
predicted.  A prediction  of  marginal  to  no  liquefaction  is  in  between  the  predictions  obtained  by 
the  SPT-  and  Vs-based  procedures  (see  Figs.  3.21,  3.22  and  3.23). 
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Corrected  CPT  Tip  Resistance,  qclN 


Fig.  3.24  - Comparison  of  Liquefaction  Assessment  Chart  Based  on  qclN  and  CSR 

Recommended  by  the  NCEER  Workshop  (Youd  et  al.,  1997;  after  Robertson  and 
Wride,  1997)  with  Results  from  the  Approach  to  Pier  Site  Between  Depths  of  6 m 
and  12  m. 
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CHAPTER  4 


SUMMARY  AND  RECOMMENDATIONS 


4.1  SUMMARY 

SASW  tests  were  conducted  along  a 240  m-long  alignment  at  the  Treasure  Island 
Approach  to  Pier  site  using  a seismic  vibrator  as  principal  source  with  receiver  spacings 
ranging  from  7.6  m to  76  m.  Additional  tests  were  conducted  at  smaller  spacings  to 
characterize  near-surface  layers  using  hand-held  hammers  and  dropped  weights  as  the  sources. 
Shear  wave  velocity  profiles  were  determined  from  the  SASW  field  data  by  forward  modeling 
using  a computer  model  with  three-dimensional  wave  propagation. 

The  zone  of  densified  sand  adjacent  to  Pier  1 at  Treasure  Island  was  correctly  identified 
in  VR-curves  and  Vs-profiles  obtained  from  SASW  tests.  Shear  wave  velocities  determined  for 
the  improved  soil  zone  averaged  about  25  m/s  greater  than  shear  wave  velocities  determined  for 
the  unimproved  zones.  This  represents  a 15%  increase  in  shear  wave  velocity  due  to  soil 
improvement.  Two  liquefaction  assessment  procedures  based  on  Vs  correctly  predicted  no 
liquefaction  for  the  improved  area,  and  marginal  liquefaction  for  the  unimproved  area. 
Although  less  conservative,  liquefaction  assessment  procedures  based  on  the  SPT  and  CPT 
provide  similar  predictions. 

This  study  further  supports  the  usefulness  of  in  situ  Vs  for  predicting  liquefaction 
potential,  and  demonstrates  the  potential  of  the  SASW  test  method  for  rapid  delineation  of 
weak  soil  layers.  For  large  study  areas,  a cost-effective  investigation  program  might  be  to  first 
develop  curves,  or  profiles,  of  VR  in  the  field  (assuming  an  approximate  sampling  depth  equal 
to  ^r/3  to  A,r/2).  The  VR-profiles  would  then  be  used  to  select  locations  for  determining  Vs- 
profiles,  which  require  a large  computational  effort,  and  sites  for  borehole  sampling  and 
penetration  testing. 


4.2  RECOMMENDATIONS  FOR  FUTURE  STUDY 

The  following  recommendations  are  provided  to  identify  areas  that  need  further  study. 

1 . Develop  a better  measure  of  goodness  of  fit  between  the  experimental  and  theoretical 
dispersion  curves.  The  maximum  likelihood  method  formulation  used  in  this  study  calculated 
the  root  mean  square  error  between  dispersion  curves.  However,  this  parameter  depends  on  the 
number  and  spacing  of  averaged  experimental  dispersion  curve  points.  A change  in  the  spacing 
of  averaged  experimental  dispersion  curve  points  in  the  asphalt  layer  can  significantly  change 
the  value  of  the  root  mean  square  error. 
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2.  Develop/improve  procedures  for  quantifying  uncertainty  in  SASW  measurements. 
Current  analytical  procedures  used  to  estimate  uncertainty  are  dependent  upon  the  variables  of 
the  starting  model.  It  would  be  useful  to  know  what  the  absolute  standard  deviation,  or 
variance,  is  for  the  Vs-profiles. 

3.  Determine  the  minimal  layer  thickness  that  can  be  accurately  resolved  at  depth.  This 
study  showed  that  thick  layers  with  different  stiffnesses  can  be  delineated  by  the  SASW  test. 
However,  it  would  be  useful  to  know  what  the  minimal  layer  thickness  that  can  be  accurately 
resolved  with  an  uncertainty  in  Vs  of  ±5%  at  a depth  of  say  10  m. 

4.  Compare  inversion  and  forward  modeling  computer  models  developed  by  various 
researchers  for  determining  Vs-profiles  from  surface  wave  measurements.  These  Vs-profiles 
should  then  be  compared  with  profiles  measured  by  other  seismic  methods. 

5.  Conduct  crosshole  and  other  seismic  tests  at  the  Approach  to  Pier  site  to  compare 
with  shear  wave  velocity  profiles  determined  from  SASW  testing. 

6.  Develop  general  guidelines  for  estimating  uncertainty  in  SASW  Vs-profiles. 

7.  Compile  additional  liquefaction  and  no  liquefaction  case  histories  where  SASW  tests 
have  been  performed. 


48 


APPENDIX  A 


SHEAR  WAVE  VELOCITY  PROFILES  BASED  ON  SASW  TESTS 
AT  APPROACH  TO  PIER  SITE 
ON  TREASURE  ISLAND,  CALIFORNIA 
March  23-24, 1996 
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Table  A.l  - Shear  Wave  Velocity  Profile  Data  from  SASW  Tests  at  Approach  to  Pier  Site. 


Assumed  Values 

Layer 

No. 

Layer 

Thickness 

Layer  Depth 

P-Wave 

Velocity 

Poisson's 

Ratio 

Density 

Damping 

Estimated 

S-Wave 

(m) 

(m) 

(m/s) 

(Mg/m3) 

Velocity  a 

(m/s) 

(a) 

Test  Array  SA-1  (unimproved  area;  STA  0+00  m to  STA  0+30  m) 

l 

0.16 

0.00  to  0.16 

— 

0.20 

1.92 

0.02 

1615 

2 

0.3 

0.16  to  0.5 

— 

0.33 

1.76 

0.02 

152 

3 

0.6 

0.5  to  1.1 

— 

0.33 

1.76 

0.02 

131 

4 

0.9 

1.1  to  2.0 

— 

0.33 

1.76 

0.02 

143 

5 

2.3 

2.0  to  4.3 

1524 

— 

1.92 

0.02 

155 

6 

3.9 

4.3  to  8.2 

1524 

— 

1.92 

0.02 

162 

7 

4.6 

8.2  to  12.8 

1524 

— 

1.92 

0.02 

174 

8 

3.0 

12.8  to  15.8 

1523 

— 

1.92 

0.02 

180 

9 

11.6 

15.8  to  27.4 

1524 

— 

1.92 

0.02 

180 

10 

15.2 

27.4  to  42.7 

1524 

— 

1.92 

0.02 

168 

11 

42.7 

42.7  to  85.3 

1524 

— 

1.92 

0.02 

274 

12 

— 

Half-space 

— 

0.33 

1.92 

0.02 

1829b 

(b)  Test  Array  SA-2  (unimproved  area;  STA  0+30  m to  STA  0+60  m) 


1 

0.16 

0.00  to  0.16 

— 

0.20 

1.92 

0.02 

1615 

2 

0.3 

0.16  to  0.5 

— 

0.33 

1.76 

0.02 

152 

3 

0.6 

0.5  to  1.1 

— 

0.33 

1.76 

0.02 

131 

4 

0.9 

1.1  to  2.0 

— 

0.33 

1.76 

0.02 

143 

5 

2.3 

2.0  to  4.3 

1524 

— 

1.92 

0.02 

155 

6 

3.9 

4.3  to  8.2 

1524 

— 

1.92 

0.02 

162 

7 

4.6 

8.2  to  12.8 

1524 

— 

1.92 

0.02 

174 

8 

3.0 

12.8  to  15.8 

1523 

— 

1.92 

0.02 

183 

9 

11.6 

15.8  to  27.4 

1524 

— 

1.92 

0.02 

180 

10 

15.2 

27.4  to  42.7 

1524 

— 

1.92 

0.02 

168 

11 

42.7 

42.7  to  85.3 

1524 

— 

1.92 

0.02 

274 

12 

— 

Half-space 

— 

0.33 

1.92 

0.02 

1829b 

(c)  Test  Array  SA-3  (unimproved  area;  STA  0+60  m to  STA  0+90  m) 


1 

0.16 

0.00  to 

0.16 



0.20 

1.92 

0.02 

1615 

2 

0.3 

0.16  to 

0.5 

— 

0.33 

1.76 

0.02 

155 

3 

0.6 

0.5  to 

1.1 

— 

0.33 

1.76 

0.02 

128 

4 

0.9 

1.1  to 

2.0 

— 

0.33 

1.76 

0.02 

137 

5 

2.3 

2.0  to 

4.3 

1524 

— 

1.92 

0.02 

140 

6 

3.9 

4.3  to 

8.2 

1524 

— 

1.92 

0.02 

168 

7 

4.6 

8.2  to 

12.8 

1524 

— 

1.92 

0.02 

177 

8 

3.0 

12.8  to 

15.8 

1523 

— 

1.92 

0.02 

186 

9 

11.6 

15.8  to 

27.4 

1524 

— 

1.92 

0.02 

189 

10 

15.2 

27.4  to 

42.7 

1524 

— 

1.92 

0.02 

162 

11 

42.7 

42.7  to 

85.3 

1524 

— 

1.92 

0.02 

271 

12 

— 

Half-space 

— 

0.33 

1.92 

0.02 

1829b 

aBased  on  3-D  computer  model  described  by  Roesset  et  al.  (1991)  and  maximum  likelihood  method 
formulation  by  Joh  (1996). 

b Assumed  value,  based  on  shear  wave  velocities  measured  in  the  Franciscan  bedrock  at  the  fire  station  and 
on  Yerba  Buena  Island  (de  Alba  and  Faris,  1996b). 
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Table  A.l  (cont.)  - Shear  Wave  Velocity  Profile  Data  from  SASW  Tests  at  Approach  to  Pier  Site. 


Assumed  Values 

Layer 

No. 

Layer 

Thickness 

Layer  Depth 

P-Wave 

Velocity 

Poisson's 

Ratio 

Density 

Damping 

Estimated 

S-Wave 

(m) 

(m) 

(m/s) 

(Mg/m3) 

Velocity3 

(m/s) 

(d)  Test  Array  SA-4a  (unimproved  area;  STA  0+90  m to  STA  1+20  m) 


1 

0.16 

0.00  to  0.16 

- 

0.20 

1.92 

0.02 

1615 

2 

0.3 

0.16  to  0.5 

- — 

0.33 

1.76 

0.02 

149 

3 

0.6 

0.5  to  1.1 

— 

0.33 

1.76 

0.02 

128 

4 

0.9 

1.1  to  2.0 

— 

0.33 

1.76 

0.02 

137 

5 

2.3 

2.0  to  4.3 

1524 

- — 

1.92 

0.02 

143 

6 

3.9 

4.3  to  8.2 

1524 

— 

1.92 

0.02 

174 

7 

4.6 

8.2  to  12.8 

1524 

- — 

1.92 

0.02 

177 

8 

3.0 

12.8  to  15.8 

1523 

— 

1.92 

0.02 

189 

9 

11.6 

15.8  to  27.4 

1524 

• — 

1.92 

0.02 

186 

10 

15.2 

27.4  to  42.7 

1524 

— 

1.92 

0.02 

165 

11 

42.7 

42.7  to  85.3 

1524 

— 

1.92 

0.02 

271 

12 

— 

Half-space 

— 

0.33 

1.92 

0.02 

1829b 

(e) 

Test  Array  SA-4b  (unimproved  area;  STA  1+03  m to  STA  1+33  m) 

1 

0.16 

0.00  to  0.16 

— 

0.20 

1.92 

0.02 

1615 

2 

0.3 

0.16  to  0.5 

— 

0.33 

1.76 

0.02 

155 

3 

0.6 

0.5  to  1.1 

— 

0.33 

1.76 

0.02 

128 

4 

0.9 

1.1  to  2.0 

— 

0.33 

1.76 

0.02 

134 

5 

1.1 

2.0  to  3.0 

1524 

— 

1.92 

0.02 

152 

6 

1.2 

3.0  to  4.3 

1524 

— 

1.92 

0.02 

146 

7 

1.5 

4.3  to  5.8 

1524 

— 

1.92 

0.02 

168 

8 

2.4 

5.8  to  8.2 

1523 

— 

1.92 

0.02 

168 

9 

4.6 

8.2  to  12.8 

1524 

_ — 

1.92 

0.02 

180 

10 

3.0 

12.8  to  15.8 

1524 

- — 

1.92 

0.02 

189 

11 

11.6 

15.8  to  27.4 

1524 

— 

1.92 

0.02 

186 

12 

15.2 

27.4  to  42.7 

1524 

- — 

1.92 

0.02 

165 

13 

42.7 

42.7  to  85.3 

1524 

— 

1.92 

0.02 

271 

14 

— 

Half-space 

— 

0.33 

1.92 

0.02 

1829b 

aBased  on  3-D  computer  model  described  by  Roesset  et  al.  (1991)  and  maximum  likelihood  method 
formulation  by  Joh  (1996). 

b Assumed  value,  based  on  shear  wave  velocities  measured  in  the  Franciscan  bedrock  at  the  fire  station  and 
on  Yerba  Buena  Island  (de  Alba  and  Faris,  1996b). 
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Table  A.  1 (cont.)  - Shear  Wave  Velocity  Profile  Data  from  SASW  Tests  at  Approach  to  Pier  Site. 


Assumed  Values 

Layer 

No. 

Layer 

Thickness 

Layer  Depth 

P-Wave 

Velocity 

Poisson's 

Ratio 

Density 

Damping 

Estimated 

S-Wave 

(m) 

(m) 

(m/s) 

(Mg/m3) 

Velocity3 

(m/s) 

(f)  Test 

Array  SA-4c  (unimproved  area;  STA  1+03  m to  STA  1+33  m;  offset  6.1  m towards  bay,  see  Figure  3.2) 

1 

0.16 

0.00  to  0.16 

— 

0.20 

1.92 

0.02 

1615 

2 

0.3 

0.16  to  0.5 

— 

0.33 

1.76 

0.02 

149 

3 

0.6 

0.5  to  1.1 

— 

0.33 

1.76 

0.02 

119 

4 

0.9 

1.1  to  2.0 

— 

0.33 

1.76 

0.02 

128 

5 

1.1 

2.0  to  3.0 

1524 

— 

1.92 

0.02 

140 

6 

1.2 

3.0  to  4.3 

1524 

— 

1.92 

0.02 

153 

7 

1.5 

4.3  to  5.8 

1524 

— 

1.92 

0.02 

158 

8 

2.4 

5.8  to  8.2 

1523 

— 

1.92 

0.02 

162 

9 

4.6 

8.2  to  12.8 

1524 

— 

1.92 

0.02 

186 

10 

3.0 

12.8  to  15.8 

1524 

— 

1.92 

0.02 

192 

11 

11.6 

15.8  to  27.4 

1524 

— 

1.92 

0.02 

183 

12 

15.2 

27.4  to  42.7 

1524 

— 

1.92 

0.02 

165 

13 

42.7 

42.7  to  85.3 

1524 

— 

1.92 

0.02 

274 

14 

— 

Half-space 

— 

0.33 

1.92 

0.02 

1829b 

(g 

) Test  Array  SA-5b  (improved  area;  STA  1+20  m to  STA  1+50  m) 

1 

0.05 

0.00  to  0.05 

— 

0.20 

1.92 

0.02 

1439 

2 

0.04 

0.05  to  0.09 

— 

0.20 

1.92 

0.02 

945 

3 

0.06 

0.09  to  0.15 

— 

0.20 

1.92 

0.02 

512 

4 

0.3 

0.15  to  0.5 

— 

0.33 

1.76 

0.02 

244 

5 

0.6 

0.5  to  1.1 

— 

0.33 

1.76 

0.02 

241 

6 

0.9 

1.1  to  2.0 

— 

0.33 

1.76 

0.02 

195 

7 

2.3 

2.0  to  4.3 

1524 

— 

1.92 

0.02 

174 

8 

3.9 

4.3  to  8.2 

1523 

— 

1.92 

0.02 

189 

9 

4.6 

8.2  to  12.8 

1524 

— 

1.92 

0.02 

204 

10 

3.0 

12.8  to  15.8 

1524 

— 

1.92 

0.02 

192 

11 

11.6 

15.8  to  27.4 

1524 

— 

1.92 

0.02 

180 

12 

15.2 

27.4  to  42.7 

1524 

— 

1.92 

0.02 

177 

13 

42.7 

42.7  to  85.3 

1524 

— 

1.92 

0.02 

286 

14 

— 

* Half-space 

— 

0.33 

1.92 

0.02 

1829b 

aBased  on  3-D  computer  model  described  by  Roesset  et  al.  (1991)  and  maximum  likelihood  method 
formulation  by  Joh  (1996). 

b Assumed  value,  based  on  shear  wave  velocities  measured  in  the  Franciscan  bedrock  at  the  fire  station  and 
on  Yerba  Buena  Island  (de  Alba  and  Faris,  1996b). 
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Table  A.  1 (cont.)  - Shear  Wave  Velocity  Profile  Data  from  SASW  Tests  at  Approach  to  Pier  Site. 


Assumed  Values 

Layer 

No. 

Layer 

Thickness 

Layer  Depth 

P-Wave 

Velocity 

Poisson's 

Ratio 

Density 

Damping 

Estimated 

S-Wave 

(m) 

(m) 

(m/s) 

(Mg/m3) 

Velocity3 

(m/s) 

(h)  Test  Array  SA-6  (improved  area;  STA  1+50  m to  STA  1+80  m) 


1 

0.05 

0.00  to 

0.05 



0.20 

1.92 

0.02 

1439 

2 

0.04 

0.05  to 

0.09 

— 

0.20 

1.92 

0.02 

945 

3 

0.06 

0.09  to 

0.15 

— 

0.20 

1.92 

0.02 

512 

4 

0.3 

0.15  to 

0.5 

— 

0.33 

1.76 

0.02 

244 

5 

0.6 

0.5  to 

1.1 

— 

0.33 

1.76 

0.02 

222 

6 

0.9 

1.1  to 

2.0 

— 

0.33 

1.76 

0.02 

183 

7 

1.1 

2.0  to 

3.0 

1524 

— 

1.92 

0.02 

168 

8 

1.2 

3.0  to 

4.3 

1523 

— 

1.92 

0.02 

186 

9 

1.5 

4.3  to 

5.8 

1524 

— 

1.92 

0.02 

177 

10 

2.4 

5.8  to 

8.2 

1524 

- — 

1.92 

0.02 

195 

11 

4.6 

8.2  to 

12.8 

1524 

— 

1.92 

0.02 

207 

12 

3.0 

12.8  to 

15.8 

1524 

— 

1.92 

0.02 

192 

13 

11.6 

15.8  to 

27.4 

1524 

— 

1.92 

0.02 

180 

14 

15.2 

27.4  to 

42.7 

1524 

■ — 

1.92 

0.02 

177 

15 

42.7 

42.7  to 

85.3 

1524 

- — 

1.92 

0.02 

286 

16 

— 

Half-space 

— 

0.33 

1.92 

0.02 

1829b 

(i)  Test  Array  SA-7  (improved  area;  STA  1+80  m to  STA  2+10  m) 


1 

0.05 

0.00  to  0.05 

- 

0.20 

1.92 

0.02 

1439 

2 

0.04 

0.05  to  0.09 

— 

0.20 

1.92 

0.02 

945 

3 

0.06 

0.09  to  0.15 

— 

0.20 

1.92 

0.02 

512 

4 

0.3 

0.15  to  0.5 

— 

0.33 

1.76 

0.02 

247 

5 

0.6 

0.5  to  1.1 

— _ 

0.33 

1.76 

0.02 

213 

6 

0.9 

1.1  to  2.0 

- — 

0.33 

1.76 

0.02 

155 

7 

2.3 

2.0  to  4.3 

1524 

1.92 

0.02 

168 

8 

3.9 

4.3  to  8.2 

1523 

■ — 

1.92 

0.02 

186 

9 

4.6 

8.2  to  12.8 

1524 

— 

1.92 

0.02 

207 

10 

3.0 

12.8  to  15.8 

1524 

• — 

1.92 

0.02 

195 

11 

11.6 

15.8  to  27.4 

1524 

— 

1.92 

0.02 

180 

12 

15.2 

27.4  to  42.7 

1524 

— 

1.92 

0.02 

180 

13 

42.7 

42.7  to  85.3 

1524 

— 

1.92 

0.02 

283 

14 

— 

Half-space 

— 

0.33 

1.92 

0.02 

1829b 

aBased  on  3-D  computer  model  described  by  Roesset  et  al.  (1991)  and  maximum  likelihood  method 
formulation  by  Joh  (1996). 

b Assumed  value,  based  on  shear  wave  velocities  measured  in  the  Franciscan  bedrock  at  the  fire  station  and 
on  Yerba  Buena  Island  (de  Alba  and  Faris,  1996b). 
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Table  A.l  (cont.)  - Shear  Wave  Velocity  Profile  Data  from  SASW  Tests  at  Approach  to  Pier  Site. 


Assumed  Values 

Layer 

No. 

Layer 

Thickness 

Layer  Depth 

P-Wave 

Velocity 

Poisson's 

Ratio 

Density 

Damping 

Estimated 

S-Wave 

(m) 

(m) 

(m/s) 

(Mg/m3) 

Velocity2 

(m/s) 

(i) 

Test  Array  SA-8  (unimproved  area;  STA  2+10  m to  STA  2+40  m) 

1 

0.16 

0.00  to  0.16 

— 

0.20 

1.92 

0.02 

1615 

2 

0.3 

0.16  to  0.5 

— 

0.33 

1.76 

0.02 

165 

3 

0.6 

0.5  to  1.1 

— 

0.33 

1.76 

0.02 

158 

4 

0.9 

1.1  to  2.0 

— 

0.33 

1.76 

0.02 

155 

5 

2.3 

2.0  to  4.3 

1524 

— 

1.92 

0.02 

146 

6 

3.9 

4.3  to  8.2 

1524 

— 

1.92 

0.02 

162 

7 

4.6 

8.2  to  12.8 

1524 

— 

1.92 

0.02 

177 

8 

3.0 

12.8  to  15.8 

1523 

— 

1.92 

0.02 

186 

9 

11.6 

15.8  to  27.4 

1524 

— 

1.92 

0.02 

180 

10 

15.2 

27.4  to  42.7 

1524 

— 

1.92 

0.02 

177 

11 

42.7 

42.7  to  85.3 

1524 

— 

1.92 

0.02 

290 

12 

— 

Half-space 

— 

0.33 

1.92 

0.02 

1829b 

aBased  on  3-D  computer  model  described  by  Roesset  et  al.  (1991)  and  maximum  likelihood  method 
formulation  by  Joh  (1996). 

b Assumed  value,  based  on  shear  wave  velocities  measured  in  the  Franciscan  bedrock  at  the  fire  station  and 
on  Yerba  Buena  Island  (de  Alba  and  Faris,  1996b). 
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APPENDIX  B 


DISPERSION  CURVES  FOR  SASW  TESTS 
AT  APPROACH  TO  PIER  SITE 
ON  TREASURE  ISLAND,  CALIFORNIA 
March  23-24, 1996 
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'Surface  wave  velocities  greater 
than  400  m/s,  see  Fig.  B.5a. 
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Fig.  B.l  - Comparison  of  Experimental  and  Theoretical  SASW  Dispersion  Curves  for  the 
Treasure  Island  Approach  to  Pier  Site,  Test  Array  SA-1.  (Note  that  Test  Array 
SA-1  Lies  Outside  the  Improved  Area.) 
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'Surface  wave  velocities  greater 
than  400  m/s,  see  Fig.  B.5a. 
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Fig.  B.2  - Comparison  of  Experimental  and  Theoretical  SASW  Dispersion  Curves  for  the 
Treasure  Island  Approach  to  Pier  Site,  Test  Array  SA-2.  (Note  that  Test  Array 
SA-2  Lies  Outside  the  Improved  Area.) 
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Fig.  B.3  - Comparison  of  Experimental  and  Theoretical  SASW  Dispersion  Curves  for  the 
Treasure  Island  Approach  to  Pier  Site,  Test  Array  SA-3.  (Note  that  Test  Array 
SA-3  Lies  Outside  the  Improved  Area.) 
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‘Surface  wave  velocities  greater 
than  400  m/s,  see  Fig.  B.5a. 


i i i i i i i 


I J-JL-L-LU 


1 1 L-L-LU 


I I 


0.1 


10 


100 


1000 


Wavelength,  A,R,  m 


Fig.  B.4  - Comparison  of  Experimental  and  Theoretical  SASW  Dispersion  Curves  for  the 
Treasure  Island  Approach  to  Pier  Site,  Test  Array  SA-4a.  (Note  that  Test  Array 
SA-4a  Lies  Outside  the  Improved  Area.) 
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Fig.  B.5  - Comparison  of  Experimental  and  Theoretical  SASW  Dispersion  Curves  for  the 
Treasure  Island  Approach  to  Pier  Site,  Test  Array  SA-4b.  (Note  that  Test 
Array  SA-4b  Lies  Outside  the  Improved  Area.) 
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Fig.  B.5  (cont.)  - Comparison  of  Experimental  and  Theoretical  SASW  Dispersion  Curves 
for  the  Treasure  Island  Approach  to  Pier  Site,  Test  Array  SA-4b. 

(Note  that  Test  Array  SA-4b  Lies  Outside  the  Improved  Area.) 
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’Surface  wave  velocities  greater 
than  400  m/s,  see  Fig.  B.5a. 
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Fig.  B.6  - Comparison  of  Experimental  and  Theoretical  SASW  Dispersion  Curves  for  the 
Treasure  Island  Approach  to  Pier  Site,  Test  Array  SA-4c.  (Note  that  Test 
Array  SA-4c  Lies  Outside  the  Improved  Area.) 
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‘Surface  wave  velocities  greater 
than  400  m/s,  see  Fig.  B.8a. 
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Fig.  B.7  - Comparison  of  Experimental  and  Theoretical  SASW  Dispersion  Curves  for  the 
Treasure  Island  Approach  to  Pier  Site,  Test  Array  SA-5b.  (Note  that  Test 
Array  SA-5b  Lies  Within  the  Improved  Area.) 
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Fig.  B.8  - Comparison  of  Experimental  and  Theoretical  SASW  Dispersion  Curves  for  the 
Treasure  Island  Approach  to  Pier  Site,  Test  Array  SA-6.  (Note  that  Test 
Array  SA-6  Lies  Within  the  Improved  Area.) 
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‘Surface  wave  velocities  greater 
than  400  m/s,  see  Fig.  B,8a. 
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Fig.  B.8  (cont.)  - Comparison  of  Experimental  and  Theoretical  SASW  Dispersion  Curves  for 
the  Treasure  Island  Approach  to  Pier  Site,  Test  Array  SA-6.  (Note  that  Test 
Array  SA-6  Lies  Within  the  Improved  Area.) 


68 


Surface  Wave  Velocity,  VR,  m/s 


400 


300 


200 


100 


j ^ , 

m f^r  -• 


Dispersion  Spacing 
Curve  (m) 

Cy 

Experimental,  6b 

0.15* 

\ j 

Experimental,  6b 

0.3 

A 

Experimental,  6b 

0.6 

V 

Experimental,  6b 

1.5 

O 

Experimental,  6b 

3.0 

O 

Experimental,  7 

7.6 

o 

Experimental,  7 

15.2 

A 

Experimental,  7 

30.5 

V 

Experimental,  6-7 

61.0 

o 

Theoretical 

0.15* 

□ 

Theoretical 

0.3* 

A 

Theoretical 

0.6 

V 

Theoretical 

1.5 

o 

Theoretical 

3.0 

o 

Theoretical 

7.6 

□ 

Theoretical 

15.2 

A 

Theoretical 

30.5 

V 

Theoretical 

61.0 

J— 1— L 


■1..  I , L-LL-L 


*Surface  wave  velocities  greater 
than  400  m/s,  see  Fig.  B.8a. 
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Fig.  B.9  - Comparison  of  Experimental  and  Theoretical  SASW  Dispersion  Curves  for  the 
Treasure  Island  Approach  to  Pier  Site,  Test  Array  SA-7.  (Note  that  Test  Array 
SA-7  Lies  Within  the  Improved  Area.) 
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Fig.  B.10  - Comparison  of  Experimental  and  Theoretical  SASW  Dispersion  Curves  for 
the  Treasure  Island  Approach  to  Pier  Site,  Test  Array  SA-8.  (Note  that 
Test  Array  SA-8  Lies  Outside  the  Improved  Area.) 
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